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Abstract 

The Lie- Trotter formula e A+ ^ = lirrijv^oo (e^/ N e^^ N ^j is of great utility 
in a variety of quantum problems ranging from the theory of path integrals and 
Monte Carlo methods in theoretical chemistry, to many-body and thermostatisti- 
cal calculations. We generalize it for the (/-exponential function eJx) = [1 + (1 — 
g ) x ] (i/(i-g)) ( with = e x^ and prove e<z ^4 + b + (1 - q )[AB + BA]/2) = 

N 



lim 



N^oo 



e l-(l-q)N 



(A/n) 



e l-{l-q)N 



b/n)]} 



This extended formula is ex- 



pected to be similarly useful in the nonextensive situations. 
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The Lie- Trotter product formula states that if A and B are not necessarily commuting 
finite square matrices or bounded operators with respect to some convenient norm, then 



lim 

N^oo 



(e A ' N e*l N ) 



N 



(1) 



This formula has been central in the development of path integral approaches to quantum 
theory, stochastic theory, and in quantum statistical mechanics 0,0. In particular, 
Suzuki has employed this to develop quantum statistical Monte Carlo methods || Q, 
general theory of path integrals with application to many-body theories and statistical 
physics H, and more recently to mathematical physics f|. In the past decade, the use of 
monomial form of the exponential function defined by (hereafter called g-exponential) 

q)x] 1/{1 - q) 



[1 + (1 



(qeTZ) 



(2) 



with a cut-off, for q < 1, when [1 + (1 — q)x] < (analogously, for q > 1, e q (x) diverges at 
x = l/(g— 1)), has played an important role in the development of nonextensive statistical 
physics ||. Here we consider the corresponding g-exponential operator defined similarly 
by 

e,(i) = [l + (l-g)i] 1 /(^) 



(qeTZ) 



(3) 



Here, A can be a c-number, a square matrix, or an operator. For g going to 1, this 
goes to the usual exponential form defined in Eq. (1). There is a large class of prob- 
lems (Levy-like H and correlated-like [fT(i|l anomalous diffusion, turbulence |TT[ 



m nonneu- 



tral plasma, nonlinear dynamics ||12|| , solar neutrino problem [ 13 1, cosmology [O], electron- 
positron collisions reassociation of heme-ligands in folded proteins[16], among many 
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others) with a (multi) fractal structure in the relevant space-time in which this and related 
functional forms emerge naturally. Many of these have received until now only a classical 
approach; however, their quantum counterparts could be analyzed as well, and then the 
operator form we are addressing would be the natural one to employ In this paper, we 
wish to present the generalization of the formula in Eq. (1) for these monomial functions 
in the form: 



e q (A + B + (l-q){A,B}/2) = Km { 



e l-(l-q)N 



(A/n) 



e l-(l-g)iV 



(B/N)]} N 
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where {A,B} = (Ab + BA). Note that e x ^ q)N {A/N) = [1 + (1 - q)A}^ {l ^ N \ has 

the properties that it goes to e^ A ^ N ^ both when q = 1 for finite N and for q ^ 1, but 
N — > oo, so that we recover Eq. (1) appropriately. To establish Eq. (4), we first consider 
the following pair of operators: 



C 

b 



'A + B+ {l-q){A,B}/2N\ 



N 



e q (A/N)e q (B/N) . 



(5) 



Then, following similar steps as in the case of the exponential operators in the 
norms of these operators are found to be bounded by the following inequalities: 



l|C|| < 
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and 



\D\\ < ||e 3 -(i/iV)|||| e? -(S/iV)|| 
< e 9 -(||i||/iV)e 9 -(||5||/iV) 




1 (1/(1-9)) 
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where we have used that (i) the norm of an increasing function of an operator cannot 
exceed the function of the norm of that operator, (ii) the norm of a sum of two operators 
cannot exceed the sum of the norms of those operators, and (iii) the norm of the product 
of two operators cannot exceed the product of the norms of those operators. These 
expressions are valid for bounded operators as long as the restrictions mentioned after 
Eq. (2) are obeyed. In the last line in Eq. (7), we used the identity which follows from 
Eq. (2) for c-numbers, that 



e q (x)e q (y) 



e q {x + y 



q)xy) 
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in which we set x = \\A\\/N, and y = \\B\\/N. Thus the norms of the two operators are 
found to be bounded by the same quantity. We now calculate the norm \\C N — D N \\. We 
therefore consider the operator identity given in |], || , 



N 



C N -D N = J2 C k -\C - D)D 

k=l 



N~k 



(9) 



and find a bound on the left side of this equation by the norm on the right hand side, 
using Eqs. (6) and (7): 
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Expanding C and D in a power series in (1/N), we now estimate the norm of their 
difference for large N 
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where [A, B\ = AB — BA. It is worth noticing that there is no 1/N term. 

By using inequality (11) into (10) we have thus obtained the following inequality for any 

q and large N: 
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We now let g = 1 — (1 — q)N, which implies o(-^f) = o(l/N 2 ), and take the limit N — > oo, 
thus showing that the right side of Eq. (12) tends to zero, which implies the required 
generalization of the Lie- Trotter formula, Eq. (4), i.e., 



{1 + (1 - q)[A + B + (1 - q){A, B}/2]}~ 



Urn 



N^oo 



{[1 + (1 - q)A]V=3F[l + (1 - q)B]~ 



In conclusion, we have generalized here (Eq. (4) or, equivalently, Eq. (13)) the Lie- 
Trotter formula, for g-exponential operators which occur in the development of nonex- 
tensive statistical physics |§. It is expected that this formula will be useful in quantum 
versions of generalized simulated annealing ||T^| and in accelerating the rate of conver- 
gence of algorithms used in quantum chemistry[|J, as well as for performing calculations 
of quantum nonextensive systems. 
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